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\^ [ In scenarios with a gravitino LSP, there exist strong BBN constraints on the 

abundance of a possible stau NLSP. We find that in settings with substantial left- 
• right mixing of the stau mass eigenstates these constraints can be evaded even for 

very long-lived staus. 
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1 Introduction 



Arguably, supersymmetry (SUSY) is one of the most plausible extensions of the standard 
model (SM). Apart from its theoretical appeal, SUSY has the virtue of providing a 
compelling dark matter candidate, the lightest supersymmetric particle (LSP), which is 
stable if R parity is conserved. A particular interesting dark matter candidate is the 
gravitino, which evades direct detection because all its interactions are suppressed by 
the Planck scale. The hypothesis of a gravitino LSP may nevertheless be tested at the 
LHC if certain conditions are met. First, the next-to- lightest supersymmetric particle 
(NLSP) has to be charged, and second the gravitino mass may not be too small. In this 
case, one could observe long-lived charged particles in whose decays one could probe 
the properties of the LSP [1-4]. The collider phenomenology of such scenarios has been 
explored in various studies [1-9]. 

There are several theoretical reasons which make it appear desirable to have gravitino 
masses 7713/2 not much smaller than the masses of the SM superpartners. For instance, 
simple explanations of the and B fi terms seem to require not too small 777,3/2 [10, 11]. 
Further, many simple mechanisms of baryogenesis, in particular leptogenesis [12], need 
rather high reheating temperatures Tr [13] which can be achieved for gravitino masses of 
about 10 . . . 100 GeV [14, 15] (i.e. at least the constraints from gravitino overproduction 
can be satisfied). 

However, there are severe constraints on such scenarios coming from cosmology. The 
observed primordial abundances of light elements produced in big bang nucleosynthesis 
(BBN) allow to place stringent constraints on the number density of long-lived particles 
whose decays happen during or after BBN and induce nuclear reactions that change 
the element abundances [16-18]. While a neutralino NLSP is strongly disfavored for 
gravitino masses in the GeV range [3,19], scenarios with a sneutrino NLSP are essentially 
unconstrained but very hard to test experimentally [20,21]. This makes a charged 
slepton, specifically a stau, particularly appealing as an NLSP candidate. The stau 
NLSP abundance and lifetime can satisfy the limits obtained from BBN by considering 
NLSP decays alone [22,3,15,19,23-28]. However, as pointed out in [29], charged NLSPs 
form bound states with light nuclei, which leads to a drastic overproduction of ^Li. This 
process, known as Catalyzed BBN (CBBN), leads to strong constraints on the stau relic 
abundance, unless the NLSP lifetime is shorter than a few thousand seconds. 

Several ways to circumvent BBN constraints have been discussed in the literature. 
For instance, entropy production between NLSP freeze-out and the start of BBN can 
dilute the NLSP abundance sufficiently to satisfy all constraints even for long lifetimes 
[15,30,31]. However, in order to arrive at such scenarios one usually relies on new sectors 
which are typically very hard to access experimentally. Alternatively, the NLSP can be 
sufficiently short lived if the gravitino is very light, R parity is slightly broken [32,33] or 
the superpartner mass spectrum is sufficiently heavy [34]. However in these cases it is 
practically impossible to test the nature of the LSP. 

The purpose of this study is to point out that there are regions within the parameter 
space of the minimal supersymmetric extension of the standard model (MSSM) where 
the relic stau abundance is strongly suppressed such that the bounds from CBBN can 
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be evaded even for long stau lifetimes. As we shall see, small thermal relic abundances 
of staus occur in parameter regions with a substantial left-right mixing of the stau mass 
eigenstates, where the annihilation into Higgs bosons is greatly enhanced. 

The paper is organized as follows. In the next section we will introduce the stau- 
Higgs coupling and calculate the Higgs channel cross section. Section [3] is devoted to a 
discussion of theoretical constraints on trilinear couplings between Higgs and r fields. In 
Section m we review the relevant BBN constraints and discuss the stau relic abundance. 
Continuing with Section [5] we introduce three scenarios within the MSSM in which a 
strong suppression of the stau relic abundance can be achieved such that all cosmological 
constraints can be evaded. Finally in Section [HI we briefly discuss the implications of our 
scenario for the LHC 

2 Annihilation into Higgs bosons 

In the early universe, superpartners are copiously produced; usually they are assumed 
to be in thermal equilibrium. As the universe cools down, they will cascade into staus, 
which we assume to be the NLSPs. Since staus are metastable, until the BBN era their 
abundance will only decrease due to annihilation. In most analyses performed so far, 
the lightest stau is assumed to be purely right-handed. Then, for its freeze-out, only 
electroweak annihilation processes have to be considered. The couplings governing the 
relevant reactions are either the electric charge e or the U(l)y-coupling gy = e/ cos^w, 
where denotes the Weinberg angle. These couplings are rather small, leading to a 
relatively large stau abundance after freeze-out [22], 

1%, > 10-'^ for m^, > 100 GeV , (1) 

where the abundance Y = n/s is defined as the ratio of number and entropy densities. 
Such a large relic stau abundance is allowed by CBBN only if the gravitino is very light 
and the stau lifetime accordingly short. If the lighter stau has a left-handed component, 
the electroweak annihilation cross section gets enhanced due to its SU(2)l couplings, but 
as CBBN bounds are very tight, the inclusion of further gauge interactions changes the 
situation only marginally. On the other hand, we shall see that in the case of substantial 
left-right mixing in the stau sector, the couplings between staus and Higgs bosons can 
get significantly enhanced, thus greatly suppressing the stau relic abundance. 

2.1 Coupling of staus to Higgs bosons 

To find the regions of parameter space where this annihilation reaction is important, 
we now turn to the Lagrangean term which describes the couplings between the light 
stau and the light Higgs. In our analytic discussion, we make a couple of simplifying 
assumptions; later, in Section [5l we will take into account all interactions and states. 

We shall assume that there is no generation mixing in the slepton sector which is 
suggested by flavor constraints. Furthermore we take /i and A'^ to be real parameters. 
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Then the relevant terms reaclll 

■^Tirih = 7777- \ sin (a + (3) [(tan^ 6'w - l) cos^ Or - 2 tan^ 6^ sin^ 6^] 

u COS a — A'^ sin a ^ , sin a 1 , _ , _ , , 

+ 7> sin 29r + 2ml \ ^ r^r^ . 2 

cosp cosp J 

For simphcity we assume that the Higgs bosons except h be relatively heavy (> 300GeV), 
which is the case for all models we are considering later. This allows us to work in the 
'decoupling limit' where the mixing parameter a can be written as a ^ /? — 7r/2 and 
therefore cos a ~ sin (3 and sin a 2^ — cos/3. The leading term of the Lagrangean which 
couples the hghtest stau to the Higgs is then given by 

■^rAh = sm2er {jj, tan (3 + A'^} hr^T^ 

sin 26r m~ h rj^rf . (3) 



2Mv ^i^i^ 



Here denotes the off-diagonal element of the 2x2 stau mass matrix (cf. Ap- 

pendix E| • 

2.2 Higgs channel cross section 

With the stau-Higgs coupling (IHl), we can now calculate the cross section for the annihi- 
lation of the light staus into light Higgs bosons, 

r+ + rf h + h . (4) 

The contributing Feynman diagrams are shown in Figure [H In our approximation, we 

n h Ti h 
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Ti h Ti h 

(a) (b) 

Figure 1: Stau annihilation into Higgs bosons, 
consider only the exchange of ri, neglecting the exchange of T2- In zeroth order of a 



^There exist different sign conventions for the A-parameter. Here we follow [35]. 
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velocity expansion, the thermally averaged annihilation cross section ((Tann v) is equal to 
the cross section aann times the relative velocity Wrei of the incoming staus. We obtain 



) 
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sin 26 r m~ 



m~ (2m~ — mj)"^ 



(5) 



Clearly, this annihilation cross section becomes important for sizable left-right mixing 
and relatively small stau masses. 

2.3 Comparison with electromagnetic cross section 

It is instructive to compare the cross section (I5l) with a typical electroweak cross sec- 
tion. For example the annihilation cross section of staus into photons is given by [22] 
(cann — 47ra^/m| . If m^^j and m/j are not too close, one has 



where /icfr = /i sin 26r. Hence, for yUcfr > 2mr^ the annihilation cross section is dominated 
by the Higgs channel. Even for an order one ratio yUcfr/2m?:^ one obtains a dramatic 
reduction of the primordial stau abundance. However, as we shall discuss next, there 
are constraints on the ratio fica/2m:f^, implying that the reduction cannot be arbitrarily 
strong. 

3 Theoretical constraints on /i 

The enhancement of the stau annihilation cross section relies on a large trilinear Higgs- 
stau coupling, which might lead to an unwanted (color and) charge breaking (CCB) 
minimum of the potential. Such minima might be acceptable if the 'physical' vacuum 
is sufficiently long-lived. In what follows, we will first analyze the tree-level potential 
and see that in the interesting regions of parameter space there is indeed an unphysical 
minimum. We will discuss tunneling to this vacuum. Then we will discuss quantum 
corrections to the potential and see that they lift the unwanted minimum and possibly 
even make it disappear. Finally, we will comment on constraints on fi arising from 
unitarity. 

3.1 CCB constraints at tree level 

The relevant field space is given by two real fields, r = ^Re(ri) and h. The correspond- 
ing scalar potential around the electroweak vacuum can be written as 





(6) 




(7) 
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where at tree level ruh ~ 90 GeV, bt ~ 17 GeV, Xh ^ Xrh ^ 0.018, ~ 0.028 and 
the trilinear coupling bhrr = —■^fJ'cS- For very negative bhrr values there is a second 
minimum of the tree-level potential. We have searched for minimal paths connecting the 
electroweak vacuum with the second, deeper minimum. It turns out that the relevant 
field direction is always very close to x = ^{t + h). The corresponding potential along 
that direction can be written as 

V = \{ml + ml) + - Vr^) + + + K,) . (8) 

In order to check whether the lifetime of the local minimum at x = exceeds the age of 
the universe we have to calculate the so called bounce action Sb along the lines of [36], 
which should satisfy Sb > 400 [37]. Using the potential (IHl) with the coefficients given 
by their tree-level values we find an upper bound on the coefficient of the trilinear term 
1^ - I ^ 33.5 GeV, which translates into 

yUcff < 630 GeV (for tan/3 = 50 and m^, = 120 GeV) . (9) 

Note that for /Xefr ^ 430 GeV the tree-level potential does not exhibit a second, deeper 
minimum at all. 

3.2 Quantum corrections to CCB constraints 

It is well known that quantum corrections can change the tree-level picture drastically. In 
order to analyze the situation properly, one has to take into account radiative corrections 
to the potential. This is a very complicated issue, and we refrain from performing a 
complete analysis here. In order to get a feeling for the impact of radiative corrections 
we include the standard stop loop correction (cf. e.g. [35, p. 245 f.]). The resulting 
effective potential is significantly steeper in the Higgs direction, such that the physical 
Higgs mass can violate the tree-level bound mh < mz- It also turns out that the cubic 
and quartic coefficients in the Higgs potential get enhanced by ~ 70 % and ~ 40 %, 
respectively. This has important implications for the bounce action: plugging the loop 
corrected values for these parameters into (IHl) we find that the metastable minimum has 
a sufficiently large lifetime for 

/ieff < 780 GeV (for tan/3 = 50 and m^, = 120 GeV) . (10) 

The second, deeper minimum exists only if fi^s ^ 500 GeV. 

There are similar effects, in particular for the potential in r direction. This issue has 
been studied in [38], where it was found that there are no charge breaking minima at 
all after quantum corrections are taken into account. Whether or not these statements 
also apply to parameter regions with large tan f3 and fi^s "will be studied elsewhere. 

3.3 Unitarity bound 

Further constraints on fi^s come from unitarity. The unitary cross section for scalar 
particles can be calculated using partial wave expansion. In the case of a non-elastic 
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process it takes the form [39] 
_ 47r(2J+l) 

where J is the angular momentum of the partial wave and pin is the three-momentum 
of one incoming particle. The dominant contribution to the stau annihilation into Higgs 
bosons comes from s-wave scattering. Therefore the perturbative unitarity constraint 
relevant to our discussion is 

I Pin I 

If the bound is not respected, this signals that the perturbative calculation of the cross 
section is no longer valid. In our case, the unitarity bound translates into a constraint 
on the /i parameter. The bound is practically independent of m/^, and reads 



^,._^<4.1TeVx^^j. (13) 

We find that the annihilation channel into Higgs pairs relative to the annihilation into 
gauge bosons can well be larger by factors of 0(10'^) without violating the unitarity 
bound ([H. 



4 Primordial staus 
4.1 BBN constraints 

Various cosmological constraints on the stau yield = 1~+ + 1~- have been explored 
in the literature. In a scenario where the LSP is very weakly coupled, the NLSP de- 
cays a considerable time after the start of BBN. The decay products of such long-lived 
particles can alter the primordial light element abundances. This leads to constraints 
on the released electromagnetic and hadronic energy [40-42]. The constraints on the 
stau relic abundance from decays depend on the stau lifetime and mass as well as on 
the electromagnetic and hadronic branching ratios. For reasonable values of the stau 
mass, the hadronic branching fraction is typically < 0(10"'^) which leads to rather weak 
constraints. Stronger constraints come from electromagnetic energy injection, especially 
at late times. Here the bounds can be as strong as m^^l- < 10"^^ GeV [43]. However, 
typically the stau NLSP abundance and lifetime can satisfy the limits obtained from 
BBN by considering NLSP decays alone [3, 15, 19, 23-28] 

In addition to injecting energetic showers into the plasma through decays, negatively 
charged particles can form bound states with light nuclei, which can lead to a drastic 
overproduction of ^Li [29] and ^Be [44]. This leads to strong constraints on the stau 
yield 1~- for lifetimes longer than a few thousand seconds. While [45] gives a rather 
conservative upper bound of Y^- < 10~^^ derived from ^Li alone, [44] takes into account 
^Li as well as ^Be leading to YV- < 10~^^ which translates into 1^ < 2 ■ 10~^^ in the 
absence of a stau anti-stau asymmetry. 
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4.2 Stau relic abundance 



The relic abundance of a (meta) stable particle can be calculated using the Boltzmann 
equation. For the stau yield 1^ , (again in the absence of a stau anti-stau asymmetry) it 
takes the form [46] 



where Mp = 2.43- lO^^GeV is the reduced Planck mass, ^~eq is the abundance in thermal 
equilibrium and g^. ^ 85 represents the effective number of degrees of freedom at the 
stau freeze out (cf. [46]). 

It is well known that the relic abundance of a (meta) stable particle is inversely 
proportional to its annihilation cross section times its mass (see e.g. [35]). In the case 
where the stau freeze-out is dominated by the Higgs channel, we can write the solution 
to ([11]) as 

= 10-- (^2^^) i"^) . (15) 

If stau and Higgs are not mass degenerate (m^^^ — > 5 GeV), the annihilation cross 
section ([5]) is practically independent of rrth and depends only on the stau mass m^-^ and 
the stau-Higgs coupling oc /x tan /3 sin 26r = fj.es tan /?. Our result for the yield can then 
be written as 

-15 / \5 ^1 TeV\'' / 50 



Here we neglected subleading effects of the order 10 % like e.g. Sommerfeld enhance- 
ment H Combining our result ffTUl) with the unitarity bound ffT^ we find that the lowest 
allowed abundance is y™™ ~ 4- 10"^^ for m^^^ = 120 GeV. As explained in Section [31 the 
minimal, theoretically viable abundance might well turn out to be larger than this value. 
From the conservative bound ffTOl) we infer however that nevertheless y™™ < 10~^^. 



4.3 A comment on stau asymmetries 

So far we have assumed that there is no asymmetry in the stau sector, that is, there are 
as many r"*" as r~ degrees of freedom. On the other hand, a large class of baryogenesis 
mechanisms rely on {B + L) violation by sphalerons [48,49], which leads to an excess 
of baryons over anti-baryons if there are more anti-leptons than leptons. In particular, 
leptogenesis [12] falls into this class. From this point of view it appears natural to assume 
that there is also an asymmetry in the slepton sector at the time of stau annihilation 
and freeze out. Now it is important to distinguish between slepton number conserving 
and slepton number violating annihilation processes (Figure [2]). 

^See e.g. [47] for an explanation of the Sommerfeld effect and an estimate of the errors in the general 
case of a charged relic. 
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Figure 2: Examples for slepton number conserving (a) and violating (b) annihilation processes. 



As we have seen, in settings with substantial left-right mixing in the stau sector, 
the slepton number conserving processes dominate over the violating ones. On the 
other hand, the slepton number violating processes are still effective until the stau relic 
abundance has reached a value Yr ~ 10~^^ . . . 10~^^ [22]. It is then clear that if the 
slepton asymmetry in the stau sector is of the order of the baryon asymmetry, 77^ ~ 10"^'', 
asymmetries will play no role. However, if there are order one asymmetries in the 
(s)tau sector, this might have important implications for the relic abundance of r~: in 
a situation in which there is a large excess of r"^ over r~ , the remaining r~ are more 
likely to find an annihilation partner, and hence their relic abundance can get suppressed 
more strongly. For stau lifetimes where electromagnetic bounds are not overly restrictive, 
Tr < 10'' s, a relatively large abundance of r"*", Yf:+ < 10~^^, can still be consistent with 
BBN because they do not form bound states with nuclei. 



5 Scenarios with low Yf- 

In this section we present some 'top-down motivated settings' in which the previously 
discussed strong suppression of the stau relic abundance occurs. If a scenario can have 
a large stau annihilation cross section is fully determined by the stau spectrum. The 
necessary conditions are: 

Condition 1: Substantial left-right mixing of the stau mass eigenstates Ti and T2 driven 
by a large off-diagonal stau mass matrix element Tn~^^ ~ — /i tan (3. This 
requires large /i and tan 

Condition 2: Moderate m^^^, preferably m^^j < 200 GeV. 

We will check these conditions by looking at the entries of the stau mass matrix fn~^, 
m~ and m~ as well as on the mass eigenvalues m~ and m?r„ . 

It is clear that we could simply write down the desired stau mass matrices. The 
purpose of this section however is to present soft mass patterns defined at the unification 
scale Mgut that lead to mass matrices with the above properties. To this end, it is useful 
to recall the (one-loop) renormalization group equations (RGEs) for the stau soft masses 
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(see e.g. [35]) 



d?7i~ 1 

Sn'^^ = ylS,-3gl\M,\'-g'y\M,\'--gl^SY, (17a) 

ci?Tz~ 

Stt^^^ = 2ylSr-ig'ym\' + g'ySy , (17b) 



where 



Sr = ml^+ml^^^^ + ml^^^^ + \Ar\\ (18a) 



5- 



Y 



\Y.Y,ml (18b) 



with Yj denoting the hypercharffe of the scalar i. Note that m~ ^ m~ and m~ 

* ° ° TL TL soft T-R 



2 

TRsoft 



as other contributions are tiny (cf. ( |2T1) ) 



In order to obtain large mixing (condition 1) we have to demand that the right-hand 
sides (rhs) of ( |T7I) be similar, assuming coincident stau masses at the high scale. In 
addition, we need a large off-diagonal stau mass. |/^| which determines the size of 
together with tan /3 is fixed at the weak scale by the condition of correct electroweak 
symmetry breaking which reads (at tree-level) [35] 



mjj^sm^ (3 — rri^j^ cos"^ (3 M| largo tan /3 ^ 



An unsuppressed fn-^^ can typically be realized for /i ~ 1 — 2 TeV. In principle it is not 
difficult to get ^ in this range, however one should mention here that a relatively large /i 
might be associated with a significant amount of electroweak fine-tuning as can be seen 
from ( fT9i) .In addition, very large values for /i might lead to charge breaking vacua with 
unacceptably short lifetimes (cf. Section [3l). 

The second condition is already partially fulfilled if the mixing is sizable. To further 
reduce the stau masses, the rhs of ffTTj) should not be too negative. 

In what follows we present three scenarios where the desired stau mass patterns 
arise and low relic abundances through the Higgs channel can be achieved. We use mi- 
crOMEGAs 2.0.7 [50,51] to calculate the relic abundance of the stau NLSP numerically. 
The superpartner spectrum is determined by SOFTSUSY 2.0.18 [52] whereas the Higgs 
mass is calculated with the specialized tool FeynHiggs 2.5.1 [53-56]. For the top quark 
pole mass, we use the latest best-fit value of 172.6 GeV [57]. Experimental constraints 
on the parameter space arise primarily through mass limits. We employ the LEP Higgs 
bound mh > 114 GeV and rrir^ > 100 GeV [58]. Theoretical constraints, as discussed in 
Section [3l are not shown explicitly. 

5.1 CMSSM with large tan (3 

Let us start with the constrained supersymmetric standard model (CMSSM), which is 
defined through its free parameters mi/2, mo, Aq, tan/3 and sign/i. Although we will 
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tan /5 /ni/2 [GeV] 



Figure 3: Stau relic abundance in the CMSSM. The left panel shows the dependence of the left- and 
right-handed stau masses, the off-diagonal mass and the mass eigenvalues on tan/3. In the right panel 
we plot the logarithm of the stau rehc abundance logio(^T) in the (mi/2 — tan/3) -plane. We find the 
minimal yield to be around ~ 10^^^. The white region in the right panel is excluded due to direct 
searches. 

see that the annihilation can be more efficient in other scenarios, we find it nevertheless 
worthwhile to point out that also in this scheme a major suppression of the relic abun- 
dance is possible. The important quantity here is tan j3 which we plot against the stau 
spectrum in the left panel of Figure [3] for a typical stau NLSP parameter point. 

The plot shows that and m^f^ decrease strongly with tan/3 through the tau 
Yukawa term in the RGE of because Ut oc tan/?. Since is practically inde- 

pendent of tan /3 the off-diagonal mas^ "^tlr grows like A/tan (3 which leads to a further 
reduction of m^^ through left right-mixing. However, in spite of a strong off-diagonal 
stau mass, mixing cannot get close to maximal, because of a large difference —m~^ . 
We conclude that condition 2 can easily be satisfied while condition 1 only partially. 
In summary, we obtain a significant enhancement of the stau annihilation cross section 
through the Higgs channel at large tan (3 which is however limited by the mixing angle. 
To illustrate the effect. Figure [3] shows the stau relic abundance in the CMSSM as a 
function of tan (3. 

5.2 Non- universal Higgs masses (NUHM) 

In the NUHM we can vary - in addition to the parameters of the CMSSM - the down- 
and up- type soft Higgs masses squared at the GUT scale, m|^^ and rn^j^ . We employ the 
GUT scale stability constraint m|^^ 2 l''^^! — ^ avoid dangerous vacua and electroweak 

•^We use the definition to?lr = I ™ the following. 
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Figure 4: Stau relic abundance in the NUHM. The left panel shows the dependence of the left- and 
right-handed stau masses, the off-diagonal mass and the mass eigenvalues on the soft Higgs mass mjj^. 
In the right panel we plot the logarithm of the stau relic abundance logLo(lV) for different rnj^^ and 
Some of the low r yield regions might be excluded by theoretical constraints (cf. Section [3]), 



•'H2 



Again, the white region in the right panel is excluded. 



symmetry breaking at the GUT scale [59]. It is instructive to investigate the additional 
effects on the stau spectrum compared to the CMSSM which arise through the variation 
of the soft Higgs masses: increasing m^^ leaves unchanged, but it reduces and 
m^^ dominantly through the Yukawa term in the soft mass RGEs fll7l) which contains 



mj^_^. However more interesting for us is the impact of fnjj^ which we show in the left 
panel of Figure HI We observe that a growth of reduces rrir^ and enlarges rrir^ . As 
does not appear in the Yukawa term of the RGEs ( ITTl) . this effect arises through the 



—m 



H2 



Sy-term. For the off-diagonal stau mass m^^, it is important to recall that 
at the weak scale, cf. (fT^ . Therefore, increasing rn^j^ leads to a suppression of 

The low yield parameter region is at relatively large negative rn?^^, where the GUT 
stability constraint can still be satisfied and again at large tan (3. Here both, the left-right 
mixing as well as the off-diagonal stau mass Jt^tlj^, can further be enhanced compared 
to the CMSSM case. Stau masses which are slightly larger can be reduced by a positive 
"^j^. We conclude that both conditions can be satisfied in this region, leading to an 



extremely suppressed stau abundance. 



5.3 Scenarios with non-universal gaugino masses (NUGM) 

The possibility of having non- universal gaugino masses as high scale boundary conditions 
even for unified theories has been realized long ago [60]. For concreteness we focus 
here on SU(5) GUTs and assume that supersymmetry be broken by F-term vacuum 
expectation values of chiral fields. These fields have to transform as the symmetric 
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[GeV] mi/2 [GeV] 

Figure 5: Stau relic abundance in the NUGM. The left panel shows the dependence of the left- and 
right-handed stau masses, the off-diagonal mass and the mass eigenvalues on TO24- In the right panel 
we plot the logarithm of the stau rehc abundance logio(^) (™i/2 — "I24) -plane. Some of the 

low T yield regions might be excluded by theoretical constraints (cf. Section Once more, the white 
region in the right panel is excluded. 

product of two adjoint representations of the GUT group, but not necessarily as singlets. 
In the following we concentrate on the two smallest possible representations for the 
supersymmetry breaking fields, which are simply the singlet and the 24plet. 

The high-scale mass patterns of the gauginos of SU(3)c, SU(2)l and U(l)y turn out 
to be given as linear combinations of singlet (mi/2) and 24plet (77124) contributions [61], 

Ml = mi/2 - 0.5 77124 , 
M2 = 77ii/2 - 1.5 77724 , 

M3 = 7711/2 + 77724. (20) 

The only new parameter compared to the CMSSM is the mass arising from the 24plet, 
77724- Its impact on the stau spectrum is illustrated in the left panel of Figure [51 

777^:^ decreases quickly for growing 7/124 as the smaller M2 reduces the SU(2)l gauge 
contribution to 777-^ increases slightly, although the U(l)y gauge contribution 

to m~ shrinks. This is because the reduction of m~ through the Yukawa term in 

J r. soft T-Rsoft ^ 

( ]17b[) is less effective for smaller 7?7p^ The off-diagonal mass matrix entry 777-^^ oc n 
grows slowly as fi gets a contribution from the increasing gluino mass M3. 

It turns out that diagonal and off-diagonal stau soft masses can come very close, 
leading to maximal left-right mixing in the stau sector with a strong reduction of 771:^^ 
driven by a large TrLr^^^. As this is the key ingredient for strongly enhanced couplings of 
Ti to the light Higgs, the stau relic abundance becomes very small, as can be seen in the 
right panel of Figure [5l 
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5.4 Further remarks 

Let us briefly summarize the main results of this section. We have seen that a strong 
suppression of the stau relic abundance can be achieved such that all cosmological con- 
straints can be evaded. We also checked that the shown regions of parameter space 
are not in conflict with precision measurements of — 57 and the Tevatron limit of 
Bs — ^ fi~^fi~. Furthermore, the parameter regions with low relic abundance are all within 
the 2 a interval of the measured anomalous magnetic moment of the muon, the only ex- 
ception being the CMSSM for mi/2 ^ 1100 GeV. Comparing the low yield regions with 
the CBBN bound < 2 ■ 10"^^ we see that even in the CMSSM there remains some 
viable region in parameter space while for the NUHM and NUGM the stau relic abun- 
dance can be even smaller. Further analysis of the experimental signatures of such low 
stau yield regions appears desirable. 

6 Prospects for the LHC 

A very appealing property of our scenario is that it will be tested at the LHC. Given 
the stau spectrum, the stau-Higgs coupling is fully determined. It is then easy to see 
the impact of the Higgs channel for stau annihilation in the early universe. 

A common feature of the low yield regions is a rather small stau mass, m:^^ < 250 GeV, 
while the other SUSY particles may be quite heavy. In general, the prospects for the 
LHC depend mainly on the mass scale of the colored superpartners. If they are not too 
heavy, as can be the case e.g. in the NUHM, they will be produced in large numbers at 
the LHC due to their large cross sections. In this case one has a good chance to measure 
the whole SUSY spectrum in the cascade decays of gluinos and squarks and it should 
not be too difficult to extract information about the stau-Higgs coupling. 

If, however, the mass scale of colored particles exceeds 2 — 2.5 TeV, their production 
will become very rare or even impossible at the LHC [62]. But even in this case, being 
rather light, staus could still be pair produced, e.g. in the Drell-Yan process qq 
through a virtual photon or Z boson [63]. The number of produced stau pairs for 
different stau masses are shown in Table [TJ Further details on the stau spectrum might 
be extracted from the Drell-Yan production of ■fi'v^^ and T2T2, if kinematically 
accessible. 



Mass 


150 GeV 


200 GeV 


300 GeV 


400 GeV 




2000 


800 


200 


60 




6000 


2000 


450 


150 



Table 1: Estimated number of produced stau pairs in Drell-Yan processes at the LHC for integrated 
luminosity of 100 fb^^ extracted from Figure 1 of [64]. The number of produced mass eigenstates 
depends on the mixing angle and should lie in between the values given for left- and right-handed staus. 

Another interesting pair production mechanism in our scenario is gluon-gluon fusion. 
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Here two gluons generate a fermion loop (preferably a top-loop) to which a virtual 
Higgs boson is attached which finally decays into a stau pair. In an early study [65] 
the cross section for this process was found to be three orders of magnitude below the 
corresponding Drell-Yan cross section. Note that in our scenario, due to the strong 
stau-Higgs coupling, this suppression can at least partially be compensated such that 
gluon-gluon fusion could be comparable to the Drell-Yan process. 

Altogether we see that even in scenarios in which some superpartners are beyond the 
reach of the LHC, one may nevertheless establish the existence of supergravity in nature 
along the lines of [1-9]. 



7 Conclusions 

We have analyzed stau NLSP scenarios. In contrast to previous studies, we have not 
assumed that the stau mass eigenstates be purely right- or left-handed, but have allowed 
for non-trivial left-right mixing. In the case of substantial mixing, the annihilation into 
Higgs bosons can dominate over other channels, such that the thermal relic stau abun- 
dance, i.e. the abundance obtained without invoking late-time entropy production, can 
be strongly reduced. This makes it possible to evade all BBN constraints. The emerging 
scenarios have the advantage that they allow for rather large reheating temperatures, 
as required for instance by leptogenesis, and the cold dark matter can be explained in 
terms of 'thermally' produced gravitinos. Most importantly, all ingredients of our low 
stau yield scenarios will be tested at the LHC. 
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A Stau Masses & Mixings 

Let us briefiy introduce our conventions concerning the masses and mixings of the stau. 
We assume that there is no mixing between different slepton generations and take fi and 
to be real parameters, such the stau mass matrix can be written as 



m~ . + (sin^^w - |)^|cos2/3 + m^ -■m^{A^ + fi tan (3) 



^ ' -m^(A^ + /itan/3) m| - sin^ ^w^f cos 2/3 + 



2 2 
2 2 

m~ m~ 



TR soft 

(21) 
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A non-zero off-diagonal element fn-^^ leads to a mixing of the chiral states tl and tr. 



We can diagonalize the stau mass matrix by an orthogonal transformation 

The mass eigenvalues are 
1 



m| +m| :^ . (ml - ml Y + Am^ . (23) 



The orthogonal 2x2 matrix O is parameterized by the stau left-right mixing angle Or-, 
which relates the mass eigenstates and the chiral states, 



(24) 



Tl \ _ / cos^^^ sin^^:\ / Tl 

T2 J \— sin Or COS Or) \ 

The mixing angle 9^ is given by 

2 

—m~ 

cose, = I . (25) 

m~ — m~ )2 + mi 
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